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Abstract 



In response to a sound stimulus, the inner ear emits sounds called otoacous- 
tic emissions. While the exact mechanism for the production of otoacoustic 
emissions is not known, active motion of individual hair cells is thought to 
play a role. Two possible sources for otoacoustic emissions, both localized 
within individual hair cells, include somatic motility and hair bundle motil- 
ity. Because physiological models of each of these systems are thought to be 
poised near a Hopf bifurcation, the dynamics of each can be described by 
the normal form for a system near a Hopf bifurcation. Here we demonstrate 
that experimental results from three-frequency suppression experiments can 
be predicted based on the response of an array of noninteracting Hopf os- 
cillators tuned at different frequencies. This supports the idea that active 
motion of individual hair cells contributes to active processing of sounds in 
the ear. Interestingly, the model suggests an explanation for differing results 
recorded in mammals and nonmammals. Key words: Hair Cells; Inner Ear; 
Hopf Bifurcation; Otoacoustic Emissions 
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Introduction 

The inner ear is more than a passive recorder of sounds. It also actively 
processes sounds using metabolic energy to spectrally analyze and amplify 

the stimulus aim. One consequence of the inner ear's active 

sound processing is that it produces sounds called otoacoustic emissions. 
Otoacoustic emissions, which consist of combinations of sounds at discrete 
frequencies, can occur either in the absence or in the presence of a sound 



stimulus 



ft 



f). The exact mechanism responsible for the active processing of 
sounds and the relatedproduction of otoacoustic emissions within the ear is 
not well known (j2j, y, |7|) . Recording the emissions spectrum provoked by a 
stimulus provides a way to probe the physiological systems responsible for 
active processing of sound. 

In nonmammals, active sound processing is thought to occur within indi- 
vidual hair cells (8|, |9|). Hair cells are mechanotransduction cells responsible 



for translating sound-induced mechanica 



motion into an electrical signal 



Each hair cell consists of a 



that is received by the auditory nerve 
cell body which is contacted by the auditory nerve and a bundle of actin- 
supported fibers called stereocilia. When sound stimulates the auditory 
organ, the hair bundle is set into motion, causing transduction channels to 
be mechanically pulled open. Potassium ions flow through the transduction 
channels depolarizing the cell and ultimately causing the firing of the au- 
ditory nerve. In nonmammals, each hair cell responds preferentially at a 
specific frequency, a quality that makes the hair cell a prime suspect in the 
search for the source of the discrete- frequency otoacoustic emissions. 
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Active motion of the hair bundle is considered to be a possible mechanism 
for active sound processing in both the mammalian and nonmammalian 



car 



ia 



13 



14, 



15 



16, 



17). Experiments have shown that the hair 



bundle responds with more energy than the stimulus energy if stimulated 



near its resonance frequency (Il8l ). It has been proposed that when the hair 



bundle is displaced, calcium enters through the transduction channels and 



binds to a site inside the hair bundle ()18 



191 ) . This binding causes a change 



in the tension of the transduction channels which results in the motion of the 
hair bundle. In mammals, there is another source of active hair cell motion. 
In response to depolarization, the cell bodies of outer hair cells contract due 
to the action of the protein prestin 

Either the hair bundle motility or the outer hair cell somatic motility 
could be involved in the production of otoacoustic emissions. Interestingly, 
a physiologically-based model for hair bundle motion has been shown to 
be poised near a Hopf bifurcation for physiologically reasonable parame- 



ters (1 14 ). The motion of the outer hair cells also displays a resonance re- 



sp „n 5 e m that „ suspected t o „ise tan a physical system that is 
tuned near a Hopf bifurcation (|24l ) . 

Assuming both the hair cell bundle motion and the outer hair cell mo- 
tion is produced by a system poised near a Hopf bifurcation, the dynamics 
either system can be described by the normal form for a system near a Hopf 



bifurcation (|25l ). 



dA 
~dt 



(a + ib)A- (c + id)\A\ 2 A. 



(1) 
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The response properties of Eq. Q]have been shown to reproduce qualitatively 



23). 



many of the amplification and tuning properties of the inner ear (26, 

The otoacoustic emissions produced by the ear in response to multifre- 
quency stimuli provide ample data concerning the active processing prop- 
erties of the inner ear J(]). Here, we consider the predictions of the Hopf 
oscillator model for three-frequency forcing experiments. It is of interest to 
determine whether observed otoacoustic emissions can be explained by an 
array of Hopf oscillators, each modeled by Eq. Q] and, if so, whether cou- 
pling between the motion of the oscillators is required to obtain observed 
otoacoustic emissions results. We find that an array of noninteracting Hopf 
oscillators, perhaps describing the motion of the hair bundles or outer hair 
cells, is adequate to qualitatively explain the results of the three-frequency 
forcing experiments in both nonmammals and mammals. 



Analysis 

Assuming both the motion of the hair bundle and the motion of the outer 
hair cell body can be modeled by a system tuned near a Hopf bifurcation, 
the dynamics of each can be described by the normal form for a system 
near a Hopf bifurcation, Eq. [TJ In the normal form, the parameter a is a 
measure of proximity to the bifurcation point. When a is small in magnitude 
and negative, the 'cell' is tuned slightly below the Hopf bifurcation and 
responds to brief disturbances with decaying oscillations. If a is greater 
than zero, the 'cell' is tuned above the Hopf bifurcation and the hair bundle 
oscillates spontaneously. The parameter b is the natural frequency of the 
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cell at the onset of oscillation and d is a measure of the shift in the frequency 
of the cell as the response amplitude increases. The parameter c determines 
whether the system is supercritical (c > 0) or subcritical (c < 0). Here, 
we will concentrate on the supercritical case because it allows for small- 
amplitude, spontaneous oscillations near the bifurcation point similar to the 



spontaneous hair bundle oscillations that are observed experimentally (|28l). 



If a small time-dependent forcing is applied the system (|27l . 
normal form must be modified to include a forcing term, F, 



29, 



30), the 



d A 

— = (a + ib)A - (c + id)\A\ 2 A + F. 



(2) 



In the case of single-frequency forcing, F = fe lujt , the system can be ana- 
lyzed by considering hair bundle motions responding at the same frequency 
as the forcing frequency. Substituting A = R e tujt +' lc l ) into Eq. yields the 
following simple relationship between forcing amplitude and response am- 
plitude, 



(aR - cR 3 ) 2 + ([& - uj]R - dR 



3^2 



(3) 



271 ) each 



Using this relationship, Egufluz et al. and Camalet et al. (|26l . 
demonstrated that a generic system poised near a Hopf bifurcation displays 
many of the amplification and tuning properties that are observed in the 
auditory system. 

Two-frequency forcing experiments have been useful in studying the 
properties of otoacoustic emissions and determining their source. In sup- 



pression experiments (|31 



320. the cochlea is stimulated by a primary tone 
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as well as a second softer tone, referred to as a suppressor tone. The addition 
of the softer tone has an effect on the magnitude of the cochlear response 
at the primary frequency. Specifically, as the frequency of the suppressor 
tone approaches the frequency of the primary tone, the magnitude of the 
component of the otoacoustic emission at the primary tone decreases. The 
biological interpretation of this is that since the maximum suppression oc- 
curs when the suppressor tone is near the primary frequency, it is likely 
that the otoacoustic emission originates near the part of the cochlea tuned 
at the primary frequency. Analysis of a Hopf oscillator tuned at the pri- 
mary frequency and forced by a primary and suppressor tone supports the 



biological interpretation. Recently, Stoop et al. (|33l ). by analyzing a Hopf 
oscillator model showed that the effect of adding a second frequency close to 
the primary frequency is to increase the effective damping of the oscillator's 
response at the primary frequency. Thus a single cell, tuned near a Hopf 
bifurcation point and near the primary frequency is adequate to reproduce 
the main qualitative features of two- frequency suppression experiments. 

When the ear is stimulated by sound containing a linear combination of 
two primary frequencies u)\ and u>2 j the otoacoustic emissions spectrum is 
more complicated to analyze because distortion product otoacoustic emis- 
sions (DPOAE's) occur at linear combinations of the stimulus frequen- 



cies da, 



34|, 



35l ). In experiments, the largest DPOAE response is observed to 
occur at the 1lo\ — L02 and 2u>2 — ^1 frequency components. The presence of 
DPOAE's allows for more complicated multifrequency forcing experiments 
in which the amplitudes of the distortion products are considered. For in- 
stance, suppression experiments can be performed in which the cochlea is 
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stimulated at a combination of two primary frequencies as well as a smaller 
amplitude suppressor tone. Then the effect of the suppressor tone on the 
response at each of the primary frequencies and the distortion product fre- 
quencies can be recorded. 

In nonmammals, multifrequency forcing experiments, including two pri- 
mary frequencies uj\ and LO2 (<^i < ^2) an d a suppressor frequency, indicate 
that maximum suppression of the 2lo\ — L02 distortion product frequency 



37 



occurs when the suppressor tone is near the uj\ frequency (|3a . 
Oddly, in mammals, the reverse trend is observed and maximum suppres- 
sion of 2wi — uj<2 occurs when the suppressor frequency is near the ui2 fre 



quency (I39l . 



40). If active hair cell motion is responsible for the production 
of otoacoustic emissions, there must be an explanation for the discrepancy 
between emissions in mammals and nonmammals. 

Here, we consider the response properties of a Hopf oscillator under 
three-frequency forcing, F = F\e Wxt + F2e tU)2t + F^e tLUit . Because the system 
is nonlinear, the response contains an infinite number of frequencies, a small 
number of which will be represented prominently. If one substitutes A = 
Aie lUlt + A2e tW2t + A-^e 1 ^ 31 into the nonlinear term from the normal form, 
|A| 2 yl, the result contains only certain frequency combinations. We will 
assume that those frequencies dominate the response, and thus consider a 
response, A, that is a linear combination of those frequency components, 
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A 



H ie iwit+i<pi _|_ giw2*+»02 _|_ ft 3e iu3t+i<t>3 



+Rii2e' 



,i(2coi —u)2)t+i4>i\2 



+ i?22ie i( ' 2u ' 2_a;i ' )t+ ^ 221 



,i(2cji-t<;3)t+ii/)ii3 , p p i(2a;2— W3)t+#223 



,i(2w 3 -a;i)t+i033i 




i(aj 2 +^3 — Wi )t+#231 



(4) 



Substituting Eq. [J] into Eq. [2] yields algebraic expressions relating the 
response amplitudes, response phases, and the forcing amplitudes, F\, F 2 , 
and F3. Frequency components not represented in Eq. H]are neglected. This 
provides an analytical description for the response of the Hopf oscillator to 
three-frequency forcing. 

Under the assumption that the cells tuned near the primary frequencies, 
u\ and ui2 and the distortion product frequencies, 2u\ — L02 and 2u>2 — 
are likely produce the greatest response at 2wi — L02, we concentrate on the 
response of those four cells. Figure [1] a.-b. shows the relationship between 
the magnitude of the component of the response at the 2loi — L02 frequency, 
-R112, and the frequency of the suppressor tone, W3, for the two cells gener- 
ated the greatest response, the cells tuned at uj\ = 300 and 2oj\ — 0J2 = 270. 
In figure Q] a.-b., as observed in the two frequency suppression case, max- 
imum suppression occurs when the suppressor frequency is tuned near the 
natural frequency of the cell. In this example, the component of the response 
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of the uj\ cell at the distortion product frequency is much louder than the 
distortion product component of the response for the other three cells. So, a 
plot of the total response of the four cells shows that maximum suppression 
occurs when the suppressor tone is tuned near u\ (figure [1] c). For larger 
values of the forcing frequency, or larger values of the nonlinear coefficients 
c and d, substantial suppression may also occur at the 2u% — 0J2 frequency 
(figure Q] d.). This result is consistent with suppression curve experiments 
in nonmammals which indicate that maximum suppression of the response 
at the distortion product frequency, 2a;i — m?. occurs when the suppressor 
frequency is near the u\ frequency 34, [3^). Some experiments also 



show a secondary dip near the distortion product frequency, as predicted by 
the model (|36l ). 

Data recorded in suppression experiments is slightly different than that 
shown in figure [I] c.-d., where the forcing amplitude was held constant for 
each curve. Typically in suppression experiments, the magnitude of forcing 
needed to reduce the component of the response at 2lo\ — u>2 by a specified 
amount is recorded as the suppressor frequency is changed. A single Hopf os- 
cillator tuned at u\ = 300 yields results similar to suppression experiments, 



34j, 



35, 



3a)- 



again with maximum suppression occurring near u\ (figure[2]) (6, 

While the Hopf oscillator model qualitatively predicts the response prop- 
erties for three-frequency suppression experiments in nonmammals, it does 
not reproduce mammalian suppression results. Recall, in mammals, it is 
observed that maximum suppression of the 2uj\ — u>2 frequency occurs when 
the suppressor tone is tuned near the L02 frequency not the uj\ frequency 
as in nonmammals. Over many trials, the Hopf oscillator model never pre- 
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dieted maximum suppression near the 0J2 frequency. The probable reason 
for the discrepancy lies in differences in physiology between mammals and 
nonmammals. In nonmammals, the hair cells are embedded in a membrane 
that lacks tuning properties, while in mammals, the hair cells are embed- 
ded in the basilar membrane (l]). The basilar membrane performs much of 
the frequency filtering in the mammalian inner ear. When sound of a given 
frequency strikes the inner ear, a traveling wave is set into motion along 
the basilar membrane. This traveling wave reaches its maximum amplitude 
at different places along the membrane depending upon the frequency of 
the stimulus. For a high frequency stimulus, the wave reaches its maximum 
amplitude closer to the base of the cochlea than it would for lower frequency 
stimulus. After the wave passes through its preferred frequency, vibrations 
at that frequency are damped. 

If the mammalian cochlea is forced at two frequencies, u\ and u>2 with 
lo\ < LO2, the hair cells tuned near the higher frequency, U2 will feel both 
frequency components of the stimuli. Because higher frequency stimuli will 
have dissipated by the time the traveling wave reaches the hair cell tuned 
at uj\, that cell will feel mainly the ui\ component of the stimulus. While 
in nonmammals, the cell tuned near lo\ is responsible for generating the 
largest portion of the distortion product otoacoustic emission, in mammals 
the cell tuned near the wi frequency does not receive the full stimulus at 
both frequency components and cannot produce as great a response at the 
distortion product frequency. Therefore, it would not be surprising if most 
of the 2uj\ — u>2 distortion product frequency was generated at the 0J2 cell 
and not the u\ cell in mammals, causing maximum suppression to occur 
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near uj2- 

Conclusions 

A model consisting of a set of noninteracting oscillators tuned near a Hopf 
bifurcation was successful in qualitatively predicting the results of three- 
frequency forcing experiments observed in mammals and nonmammals. In 
the case of nonmammals, only two Hopf oscillators tuned near uo\ and 
2lo\ — L02 were necessary to predict the results of three-tone suppression 
experiments. In mammals, a Hopf oscillator tuned near the 0J2 frequency 
correctly predicted experimental results. Which cell contributes the most is 
dictated by important differences in mammalian and nonmammalian phys- 
iology. Notably, it was not necessary to assume coupling between cells of 
different frequencies in order to qualitatively reproduce experimental data. 
Though more complicated biophysically-based models would be needed to 
produce a more quantitative agreement with the experiments, it is interest- 
ing that such a simple model can explain the main experimental features. 
These results lend support to the idea that an array of oscillators tuned 
near a Hopf bifurcation could be responsible for otoacoustic emissions and 
active sound processing in the ear. Because both the somatic motility of the 
outer hair cell and the motion of the hair bundle are thought to be modeled 
by systems poised near a Hopf bifurcation, either could play the role of the 
Hopf oscillator. 

K.A.M. is grateful for fellowship support through NSF-RTG grant 0354259. 
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Figure 1: Figures a.-b. show the magnitude of the response at the 2u>i — oo 2 
frequency, -R112, for two different cells tuned at different frequencies as the 
frequency of the suppressor tone, was varied. For a.-c. parameter values 
were set at a = -.1, c = 100, d = 100, F x = .01, F 2 = .01, F 3 = .001. 
(a) Cell 1, Natural Frequency, b = 1u x - uj 2 = 270. (b) Cell 2, Natural 
Frequency, b = u>\ = 300. (c) The total 2u)\ — uj 2 component of the response 
for four cells tuned at 270, 300, 330, and 360. (d) The total 2u% — u> 2 
component of the response for the four cells tuned at 270, 300, 330, and 
360, with c = d = 500 and other parameters the same as a.-c. 
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Figure 2: Each curve shows the amplitude of the suppressor tone, F3, needed 
to suppress the response of a single Hopf oscillator at the distortion product 
frequency, R112, by a fixed amount. For the lowest curve in the diagram, 
the response, log(i?n2) is reduced by 0.5 from its unsuppressed value. Each 
consecutive curve shows the forcing needed to reduce the response by an 
additional 0.5. Parameters were set at a = —.1, b = 300, c = 100, d = 100, 
Fi = .01, and F 2 = .01. 



